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1.  Outline  of  the  research. 

In  recent  years  the  detection  and  characterization  of  microscopic  particles 
deposited  on  a  surface  has  become  an  important  problem  both  from  the  theore¬ 
tical  and  the  experimental  point  of  view.  The  importance  of  a  practical  solution 
to  this  problem  cannot  be  overemphasized.  The  final  purpose  of  all  the  research 
that  has  been  performed  in  this  field  is,  indeed,  the  establishment  of  a  standard 
technique,  either  theoretical  or  experimental  or  both,  suitable  for  the  detection 
and  identification  of  the  particles  that  may  be  deposited  on  sufaces  that  are  sup¬ 


posed  to  be  clean. 

For  homogeneous  spherical  particles  a  good  method  of  detection  and  identi¬ 
fication  seem  to  be  the  study  of  the  characteristic  resonance  peaks  that  appear 
in  the  plot  of  the  extinction  cross  section  vs.  the  wavelength.  Nevertheless,  it 
has  been  stressed  by  several  researchers  in  the  field  that  the  extinction  spectrum 
can  imdergo  severe  modifications  when  the  particles  are  in  the  vicinity  of  an  in¬ 
terface.  This  consideration  is  the  basis  of  the  research  that  we  are  going  to  de¬ 
scribe,  because  we  were  able  to  find  a  simple  method  to  use  just  these  modifica¬ 
tions  to  gain  information  on  the  morphology  of  the  particles  of  interest. 


2.  Method  of  attack. 

The  method  that  we  used  to  perform  our  research  is  based  on  image  theory 
and  is  applicable  to  particles  in  the  vicinity  of  a  perfectly  reflecting  surface.  Let 
us  remark  that  the  assumption  of  perfect  reflectivity  is  satisfied  to  a  high  de¬ 
gree  by  almost  any  polished  metallic  surface.  Under  these  conditions  the  pro¬ 
blem  of  scattering  by  a  particle  in  the  vicinity  of  the  surface  has  been  proved  to 
be  quite  equivalent  to  the  prblem  of  the  dependent  scattering  from  the  com¬ 
pound  object  that  includes  both  the  actual  particle  and  its  image  provided  that 
the  exciting  field  is  the  superposition  of  the  field  that  comes  from  the  actual 
source  and  of  the  field  that  comes  from  the  image  source,  the  latter  field  coinci¬ 
des  with  the  field  that  is  actually  reflected  by  the  interface. 

Under  the  preceding  Contract  DAJA45-93-C-0043  we  performed  a  study  of 
the  scattering  from  anisotropic  particles  modeled  as  clusters  of  spheres  and 
were  able  to  show  that,  even  when  the  particles  are  randomly  oriented,  thei  ani¬ 
sotropy  can  be  detected  through  a  study  of  the  polarized  light  that  is  scattered 
along  the  surface.  Our  present  work  is,  in  a  sense,  a  continuation  of  the  above 


mentioned  research. 

Our  basic  idea  is  that,  if  the  extinction  spectrum  is  modified  by  the  presence 
of  a  flat  interface  this  must  be  due  to  the  boundary  conditions  at  the  interface  it¬ 
self.  Now,  the  reflection  of  a  plane  wave  has  seldom  been  studied  in  terms  of 
multipole  fields.  Therefore  we  studied  the  reflection  of  the  spherical  multipole 
fields  on  a  perfectly  reflecting  surface  and  were  able  to  show  that  the  condition 
of  reflection  imposes  certain  phase  relations  among  the  multipole  amplitudes  of 
the  incident  and  of  the  reflected  wave.  Since,  as  stated  above,  the  exciting  field 
is  just  the  superposition  of  the  incident  and  of  the  reflected  field  we  were  able 
to  show  that  several  multipole  emplitudes  of  the  exciting  field  do  vanish.  It  is 
just  this  vanishing  that  affects  the  resonance  spectrum  of  a  particle  in  the  pre¬ 
sence  of  a  perfectly  reflecting  surface. 

3.  Sketch  of  the  results. 

The  theory  on  which  ve  based  our  investigation  is  fully  explained  in  the  en¬ 
closed  paper  that  has  been  submitted  for  publication  to  the  Journal  of  the  Euro¬ 
pean  Optical  Society.  Nevertheless  it  may  be  useful  to  summarize  here  the  most 
important  results  of  our  calculatios  on  spheres  and  ag^egated  spheres. 

When  one  considers  homogeneos  spheres  it  is  possible  to  define  a  universal 
function  that  contains  all  the  information  on  the  polarization  and  on  the  direc¬ 
tion  of  incidence  of  the  exciting  field.  A  study  of  this  function  allows  one  to  pre¬ 
dict  the  behavior  of  the  resonance  peaks  for  any  spherical  scatterer  when  the 
polarization  and  the  angle  of  incidence  are  varied:  the  calculated  spectra  follow 
strictly  the  expected  behavior. 

We  also  made  some  calculations  on  the  resonance  spectra  of  aggregates  com¬ 
posed  of  two  identical  mutually  contacting  spheres.  In  this  case  the  lack  of 
spherical  symmetry  prevents  one  from  defining  any  useful  function  so  that  the 
calculations  must  be  performed  from  scratch.  A  further  difficulty  stems  from 
the  fact  that  the  multiple  scattering  processes  among  the  spheres  in  the  aggre¬ 
gate  prevant  the  resonances  to  have  a  simple  relation  to  the  resonances  of  the 
component  spheres.  For  this  reason  we  considered  only  the  resonances  from  ag¬ 
gregates  of  small  spheres  with  a  high  refractive  index.  In  this  case,  indeed,  the 
only  resonances  that  may  occur  are  those  with  1=1.  As  a  result  we  were  able  to 
perform  an  analysis  of  the  extinction  spectrum  and  to  explain  the  behavior  of 
the  observed  peaks  when  the  reflecting  surface  is  present.  In  spite  of  its  simpli¬ 
city  the  simple  model  that  we  dealt  with  gives  the  guidelines  for  the  interpreta¬ 
tion  of  the  spectra  of  more  complicated  non-spherical  particles. 
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The  effect  of  the  presence  of  a  perfectly  reflecting  surface  on  the  resonance  spectra  from  particles 
is  studied  in  the  framework  of  the  image  theory.  Due  to  the  boundary  conditions  at  the  reflecting 
surface  some  of  the  resonances  are  expected  to  disappear.  This  suppression  effect  is  studied  with 
reference  to  single  spheres  and  to  binary  aggregates  of  spheres  and  it  is  shown  that  a  careful  inves¬ 
tigation  of  the  scattered  light  as  a  function  of  the  angle  of  incidence  and  of  the  polarization  may 
give  useful  information  on  the  morphology-dependent  resonances  both  of  spherical  and  non-spherical 
particles. 


1.  Introduction 


The  electromagnetic  resonances  are  a  well  known  feature  of  the  extinction  spectrum  of  a  spherical 
scatterer.  They  occur  in  the  plot  of  the  extinction  cross  section  vs.  the  size  parameter  as  peaks 
that  both  the  Mie  theory^  and  the  experiment^”"^  permit  to  relate  to  the  multipole  amplitudes  of 
the  incident  wave.  The  analogy  between  the  radial  Helmholtz  equation  for  spheres^  and  the  radial 
Schrodinger  equation  for  spherically  symmetric  systems  led  Johnson  to  formulate  a  comprehensive 
theory  of  the  resonances  both  of  homogeneous  and  of  radially  non-homogeneous  spheres  and  to  give, 
as  a  result,  a  number  of  useful  formulas  for  the  relevant  parameters  of  each  resonance.^  Resonances 
may  also  occur  in  non  spherical  particles  but  they  cannot  be  interpreted  in  the  framework  of  the 
above  mentioned  theory  neither  for  scatterers  that  are  or  can  be  modeled  as  aggregated  spheres. 
In  the  latter  case,  indeed,  the  multiple  scattering  processes  that  occur  among  the  spheres  in  the 
aggregate  make  the  resonance  spectrum  to  be  in  no  way  a  mere  superposition  of  the  peaks  that  are 
due  to  the  component  spheres.”’^  The  resonance  spectra  proved  useful  to  gain  information  on  the 
morphology  of  the  particles^’^”^’^^  so  that  several  methods  have  been  suggested  to  simplify  their 
interpretation.  For  instance  Li  and  Chylek^^  suggested  to  modulate  the  observed  resonances  by 
illuminating  the  particle  of  interest  with  two  counterpropagating  plane  waves.  In  turn,  Johnson^^*^^ 
suggested  to  study  the  scattering  from  a  spherical  particle  coupled  to  a  perfectly  reflecting  surface. 
In  our  opinion  Johnson’s  procedure  is  rather  promising  although  the  presence  of  the  reflecting 
surface  requires  considering  the  scattering  from  a  binary  cluster  even  when  the  actual  scatterer  is 
a  single  sphere.  The  image  theory,  indeed,  states  that  the  scattering  from  a  sphere  in  the  vicinity 
of  a  perfectly  reflecting  surface  is  quite  equivalent  to  the  dependent  scattering  from  the  compound 
object  that  includes  both  the  actual  particle  and  its  image  when  illuminated  by  the  superposition 
of  the  actual  incident  field  and  of  the  field  that  comes  from  the  image  source.  Since  the  calculation 
of  the  dipendent  scattering  from  aggregated  spheres  is  a  well  established  procedure^^  whose  results 
are  in  excellent  agreement  with  the  available  experimental  data,^^d7  y^eve  able  to  use  the  image 
theory  to  produce  the  full  scattering  pattern  of  an  assembly  of  randomly  oriented  identical  clusters 
in  the  vicinity  of  a  metal  surface  for  arbitrary  direction  and  state  of  polarization  of  the  incident 
wave.^^  In  particular  we  found  that  the  analysis  of  the  polarization  of  the  scattered  light  that 
propagates  along  the  surface  may  provide  useful  information  on  the  possible  anisotropy  of  the 
particles. 

The  purpose  of  the  present  paper  is  to  show  that,  as  a  result  of  the  boundary  conditions  on  a 
perfectly  reflecting  surface,  some  of  the  multipole  amplitudes  of  the  exciting  field,  that  is  defined 
as  the  superposition  of  the  incident  and  the  reflected  field,  vanish.  As  a  consequence,  some  of  the 
expected  resonances  may  not  appear  in  the  observed  spectrum:  e.  g.  the  peaks  that  are  associated 
to  a  single  vanishing  multipole  amplitude  are  bound  to  disappear.  Using  this  suppression  effect  to 
interpretate  the  resonance  spectra  of  particles  of  arbitrary  shape  is  rather  difficult, however. 
Therefore,  our  investigation  will  deal  with  single  spheres  and  with  binary  aggregates  of  spheres 
in  the  presence  of  a  perfectly  reflecting  surface.  In  fact,  although  the  resonance  spectrum  of 
aggregated  spheres  may  be  rather  complicated,  the  guidelines  for  its  interpretation  can  be  given, 
at  least  in  simple  cases. 

In  Section  2  we  revisit  the  reflection  of  a  polarized  plane  wave  on  a  perfectly  reflecting  surface 
in  order  to  reformulate  the  problem  in  terms  of  spherical  multipole  fields. 

In  Section  3  we  discuss  the  scattering  from  particles  in  terms  of  the  transition  matrix  and  give 
a  general  definition  of  the  extinction  cross  section  that  applies  even  when  a  perfectly  refelcting 
surface  is  present. 

In  Section  4  the  resonance  spectra  from  homogeneous  spheres  and  from  aggregates  of  two  iden¬ 
tical  spheres  are  investigated  on  the  basis  of  the  theory  of  Sections  2  and  3. 

A  few  concluding  remarks  are  reported  in  Section  5. 


2,  Exciting  field 


In  the  framework  of  the  image  theory/^  the  exciting  field  is  the  superposition  of  the  actual  incident 
field  and  of  the  field  that  comes  from  the  image  source:  the  latter,  in  turn,  coincides  with  the  field 
that  is  reflected  by  the  interface  in  the  absence  of  any  scatterer.  Now,  the  reflection  of  a  plane 
wave  on  a  flat  surface  has  seldom  been  dealt  with  in  terms  of  multipole  fields. Therefore,  we 
reformulate,  in  terms  of  multipole  fields  the  relations  that  are  dictated  by  the  boundary  conditions 
between  the  amplitude  and  the  polarization  of  the  incident  and  of  the  reflected  wave.  As  a  result 
we  get  the  exciting  field  as  an  expansion  in  terms  of  spherical  multipole  fields  that  depend  on 
the  parameters  of  the  incident  wave  only.  At  first  we  deal  with  an  interface  of  general  dielectric 
properties  but  we  will  specialize  our  results  to  the  case  of  a  perfectly  reflecting  interface  when  the 
need  arises. 

We  consider  a  frame  of  reference  whose  cartesian  axes  are  individuated  by  the  unit  vectors  u^;,  Uy 
and  u^,  and  assume  that  the  halfspace  z  <  0,  the  accessible  half-space,  is  filled  by  a  homogeneous 
medium  of  refractive  index  n  while  a  different  medium  of  refractive  index  n'  fills  the  halfspace 
2:  >  0:  thus  the  interface  coincides  with  xy  plane  and  its  unit  normal  concides  with  u^.  The 
electromagnetic  plane  wave 

E/  =  E'oe/ exp[2k/  •  r],  (1) 

that  propagates  through  the  halfspace  z  <  0,  is  reflected  by  the  interface  into  the  plane  wave 

Efi  i:;  E^eRexplikR  ■  r],  (2) 

where  e/  and  gr  are  the  (unit)  polarization  vectors  of  the  incident  and  of  the  reflected  wave, 
respectively,  k/  =  nkkj  and  ki?  —  nkkR  are  the  respective  propagation  vectors,  and,  as  usual, 
k  —  cj/c.  The  time  dependence  exp(— will  be  assumed  throughout.  To  get  the  relations 
between  the  amplitudes  and  the  polarizations  of  the  incident  and  the  reflected  field  one  must  also 
consider  the  transmitted  plane  wave 

Et  —  exp[2kx  •  r]. 


with  kx  =  n'Arkx,  that  propagates  through  the  halfspace  z  >  0,  and  impose  the  boundary  condi¬ 
tions. 


{Eqgj  -f  EqBr  -j-  FqCx)  X  —  0, 

(3) 

across  the  interface.  It  is  convenient  to  introduce  two  pairs  of  unit  vectors  u/^  and  whose 

index  77  =  1,  2  distinguish  wether  they  are  parallel  (77  =  1)  or  perpendicular  (77  —  2)  to  the  plane  of 
incidence,  i.  e.  to  the  plane  that  contains  k/,  kR  and  the  z  axis.  The  orientation  is  chosen  so  that 
\iR2  =  u/2  and 

U/I  X  u/2  =  k/,  Ufii  X  Uij2  =  kfi- 

Then  Eqs.  (1)  and  (2)  can  be  rewritten  as 

E/  =  Eo  ^{ei  ■  u/,)u/,  exp[jk/  •  r] 

T] 

(4) 

Eh  E'a  ^(en  •  exp[ikH  •  r], 

(5) 

V 


and  the  application  of  the  boundary  conditions,  eq.  (3),  leads  one  to  define  the  Fresnel  coefficients 
Fjj  for  the  reflection  of  a  plane  wave  with  polarization  along  u^.  The  expression  of  the  coefficients 
Fr,  in  terms  of  the  angle  ??/  between  k/  and 
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(6a) 


F2 


cos  'dj  —  n\/ sin^  ??/ 
n'^  cos  -dj  n  \/n'^  ~  sin^  t?/ 
n  cos  ??/  —  \/n'2  —  sin^  / 
n  cos  +  y/ —  n?  sin^  t?/ 


and  their  limiting  value  for  the  case  of  a  perfectly  reflecting  surface  is 


(6b) 


F,  =  (-)^-^ 

In  terms  of  the  Fresnel  coefficients  the  relation  between  the  components  of  the  incident  and  of  the 
reflected  field  is 


Eq^ji  ‘  URrj  —  EoFj^ej  •  u/;j, 


E',  =  {-r-^EoF„ 
on  account  that 

efi  •  vlr^  -  •  u/,, 

and,  as  a  result,  the  reflected  plane  wave,  Eq.  (5),  can  be  rewritten  as 

Er  =  X!  F,,  (e/  •  Uj,  )ufiT  exp[tkfi  •  r] .  (7) 

V 

At  this  stage  we  recall  that  the  multipole  expansion  of  a  vector  plane  wave  of  wavevector  K  = 
KK  is2«-29 

E  =  i?ouexp[iK  -r]  =  EoJ]  w/^\u,K)j[^)(r,/^), 

plm 

where  we  define  the  spherical  vector  multipoles 

jW(r>  K)  =  j,{Kr)X„„  (f),  K)  =  ;^V  x  i,(A>)X,„(f), 

and  the  amplitudes 

Wl^\u,  K)  =  4«'u  •  xr„(K),  W^^\u,  K)  =  47rz'+i(K  x  u)  •  X?„(K). 

In  the  preceding  equations  the  superscripts  1  and  2  are  the  values  of  a  parity  index  p  that  distin¬ 
guishes  the  magnetic  multipoles  (p  =  1)  from  the  electric  ones  (p  =  2)  and  the  functions  'X.im  are 
vector  spherical  harmonics.^®  Using  the  notation  and  the  phase  conventions  of  Rose,^®  the  latter 
functions  can  be  defined  in  terms  of  the  spherical  harmonics  Yim  (r)  as 


X,^(f)  =  -^C(l,/. 


where 

^0  =  u^,  £±1  =  ±  juy)  (8) 

are  the  spherical  basis  vectors  and  the  symbol  C  denotes  the  Clebsch-Gordan  coefficients 
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,  x//(/  +  l)-m(mibl)  ,  m 

1)  -  C(l,/,/;0,m)  =  -^= 


Accordingly,  the  multipole  expansions  the  incident  and  the  reflected  field  are 

E;  =  Eo  J2ie,  ■  UJ,)  ^  nk)  (9a) 

Tf  plm 

E/2  =  £*0  ^  ‘  '^2 

r]  plm 

respectively,  where  we  define 

1<L  = 

In  order  to  relate  the  multipole  amplitudes  of  the  reflected  wave,  those  of  the  incident 

wave,  we  resort  to  the  spherical  basis,  Eq.  (8),  to  represent  the  vectors  as 

A  —  )^(A  •  =  y^(— 

so  that  the  dot  product  of  any  two  vectors  takes  on  the  form 
A-B  ^{-rApB.p, 


and  the  spherical  representation  of  any  unit  vector  v  whose  polar  angles  are  r?  and  (p  turns  out  to 
be28 


Now,  if  and  (pj  are  the  polar  angles  of  the  direction  of  incidence,  we  have 

ki  =  (’d],(pf),  u/i  =  (t?/ +  !,¥>/),  u/2  =  (|,y>/ +  f), 

kj{  =  (it  -  ■dj,(pj),  Ufii  =  (i9/ +  +  ?r),  Uie2  =  (f ,  y>/ +  f ), 

SO  that  it  is  an  easy  matter  to  see  that  the  properties  of  the  spherical  harmonics  under  change  of 
their  arguments  yield  the  following  relation  between  the  multipole  amplitudes  of  the  incident  and 
of  the  reflected  field 

<L-  ,  (10) 

As  a  result  the  reflected  field  takes  on  the  final  form 

Er  =  EoJ2  ■  y/-))  V(-)’'+P+'+’"  jg(r.  nArh _ (n)_ 

t]  plm 

that  depends  on  the  parameters  of  the  incident  wave  only. 

Equation  (11)  allows  us  to  write  the  exciting  field,  that  is  the  superposition  of  the  incident  and 
of  the  reflected  field,  as 

Es  =  E,  +  Eh  =  Eo  ^(e,  •  u,,)  ^[1  +  nk).  (12) 


For  a  surface  of  general  dielectric  properties  |£^|  ^  1,  so  that  the  term  within  square  brackets  in 
eq.  (12)  never  vanishes.  However,  for  a  perfectly  reflecting  interface  Fp  ~  (~)^~^*  therefore,  when 
[1  —  (— )P+^+”^]  —  0,  i.  e.  when  p-f  /  +  ^  is  even,  the  corresponding  multipole  is  not  present  in  the 
exciting  field. 
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3.  Scattered  field  and  extinction  cross  section 


The  field  that  is  scattered  by  any  particle  embedded  in  a  homogeneous  medium  of  refractive  index 
n  can  be  expanded  in  a  series  of  spherical  vector  multipoles 


where  the  multipole  fields  are  identical  to  the  multipoles  except  for  the  substitution  of 
the  spherical  Hankel  functions  of  the  first  kind,  h\^\kr)^  to  the  spherical  Bessel  functions,  ji(kr). 
The  label  rj  that  is  attached  to  F^St]  and  to  the  amplitudes  recalls  that  the  scattered  field 

depends  on  the  state  of  polarization  of  the  incident  wave.  According  to  Waterman the  multipole 
amplitudes  of  the  exciting  field  are  related  to  the  amplitudes  of  the  scattered  field  through  the 
equation 


a(p)  - 

^r)lm 


Eo(P.p')  W(P') 


where  the  quantities  are  the  elements  of  the  so  called  transition  matrix,  S,  that  accounts 

for  the  morphology  (structure  and  scattering  power)  and  the  orientation  of  the  particle.  In  the 
absence  of  any  substrate  the  amplitudes  concide  with  those  of  the  incident  plane  wave  and 

S  is  the  transition  matrix  of  the  actual  scattering  particle.  However,  when  a  perfectly  reflecting 
surface  is  present  the  amplitudes  of  the  exciting  field,  according  to  the  preceding  section,  are 

=  [1  -  (15) 

and-S  should-Te  the  transition  matrix  appr-opri-atc  to  the  compound  object  that  includes  both  the 
actual  particle  and  its  image.  In  the  latter  case,  according  to  Eq.  (15),  even  a  non  vanishing  incident 
amplitude  may  yield  a  vanishing  exciting  amplitude  and  thus  affect  the  amplitudes 

of  the  scattered  field  up  to  the  suppression  of  some  of  the  characteristic  resonance  peaks.  The 
possible  occurrence  of  this  suppression  effect  is,  in  principle,  independent  of  the  distance  of  the 
scatterer  from  the  reflecting  surface.  Nevertheless,  when  the  scatterer  that  we  consider  is  sufficiently 
far  from  the  surface  the  multiple  scattering  processes  between  the  actual  particle  and  its  image 
become  negligible  and  the  scattered  field  becomes  identical  to  the  superposition  of  the  fields  that 
are  scattered  independently  by  the  actual  particle  and  by  its  image  each  illuminated  by  the  exciting 
field. 

We  now  need  a  definition  of  the  extinction  cross  section  that  applies  even  in  the  presence  of  a 
flaF interface.  To  this  end  we  resort  to  the  optical  theorem  that,  according  to  van  de  Hulst,^  can 
be  proved  by  considering  the  field  that  is  actually  detected  by  an  optical  instrument.  Let  us  recall 
that  the  scattering  amplitude  of  any  particle,  can  be  defined  through  the  equation 


= 


exp  (inkr) 


E'of^(k5,k/), 


provided  the  particle  is  at  the  origin  and  the  distance  of  observation,  r,  is  large.  The  scattering 
amplitude  depends,  in  general,  on  the  morphology  as  well  as  on  the  orientation  of  the  scatterer 
with  respect  to  the  incident  field  and,  once  the  amplitudes  are  known,  its  expression  is 

X  X,™(ks)] .  (17) 


In  terms  of  the  optical  theorem  reads 
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47r  A  >s 

^exlr]  ”  (^5  —  ^/j^/)]) 


(18) 


where  =  fn  *  u^^,,  and  the  index  r]  recalls  that,  for  an  anisotropic  scatterer,  the  cross  section 
depends  on  the  polarization.  Now,  according  to  van  de  Hulst,^  Eq.  (18)  holds  true  even  when 
a  reflecting  surface  is  present  provided  that  the  direction  of  observation  be  the  direction  of  the 
reflected  wave  (k5  =  this  is,  indeed,  the  forward  scattering  direction  when  a  reflecting  surface 
is  present.  Nevertheless,  one  has  to  take  account  that,  according  to  Section  2,  when  the  incident 
wave  is 

E/,;  =  Eonjrjexp{iki  •  r), 

i.  e.  when  Elr,  is  either  parallel  {ij  —  1)  or  perpendicular  to  the  plane  of  incidence  (rj  =  2),  the 
reflected  wave  is 

=  Eo{-y^'^URrjexp{i'kR  •  r), 

so  that  reverses  its  phase  when  the  polarization  of  E/;^  is  perpendicular  to  the  plane  of 

incidence.  Therefore,  the  total  field  that  is  detected  by  an  optical  instrument  in  the  direction  of 
reflection  is 

Ez)r,  =  Eo[i-y~^^Rrj  exp{ikR  •r)+  f77(k/i,  k/)], 

and  the  optical  theorem  should  read 

<Te.tr,  =  (19) 

where 


Aimk  ^  ^ 


Ap)*  q(p>p')  T^(p'') 
Rrjlm  *^ltn,Vm‘  EriVm*  * 


plm  p'Vm' 


The  phase  factor  (-)'^“^  in  Eq.  (19)  does  not  appear  explicitly  in  the  expression  given  by  Johnson^^ 
because  this  author  considers  only  normal  incidence  and  includes  the  correct  phase  in  the  expression 
for  the  reflected  plane  wave. 

When  the  particle  of  interest  is  a  sphere  whose  distance  from  the  interface  is  sufficiently  large 
that  the  interaction  with  its  image  need  not  be  considered  the  exciting  field  is  still  the  superposition 
of  the  incident  and  of  the  reflected  field  but  the  transition  matrix  is  that  appropriate  to  a  single 
sphere.  In  this  case  S  is  diagonal, 

c(p.pO  -  a  (21) 


and  its  elements  are  given  by 


(1  +  n6pi)ui(nokp)u\{nkp)  -  (1  +  n6p2)u\{nokp)ui{nkp) 
(1  +  nSpi)ui{nok p)w\{nkp)  —  (1  +  nSp2)u\{nokp)wi(nkp)  ’ 


where 

n  =  — - 1,  ui(x)  —  xji{x),  wi{x)  =  xh^j^\x), 

no 

The  quantities  and  coincide  with  the  Mie  coefficients  bi  and  a/,  respectively,  for  the 
scattering  by  a  homogeneous  sphere  of  radius  p  and  refractive  index  no  embedded  in  a  homogeneous 
medium  of  refractive  index  n.  Thus,  by  considering  a  sphere  at  least  so  distant  from  the  interface, 
Eq.  (20)  can  be  rewritten  as 
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(23) 


p/m 

On  account  that  the  elements  are  independent  of  m,  it  is  meaningful  to  define  the  quantity 

=  'E^Rrl>mK1,m  (24) 

m 

whose  behavior  as  a  function  of  the  angle  of  incidence  may  give  useful  information  for  the  interpre¬ 
tation  of  the  resonance  spectrum.  In  fact  any  resonance  of  a  spherical  scatterer  that  is  associated 
to  a  vanishing  is  bound  to  disappear. 


4.  Results  and  discussion 

Before  we  go  to  discuss  the  results  of  our  specific  calculations  it  may  be  useful  to  recall  a  few 
facts  about  the  electromagnetic  resonances.  For  a  spherical  scatterer  a  resonance  occurs  when, 
with  varying  wavelength,  one  of  the  quantities  Eq.  (22),  undergoes  a  fast  change  from  a 

rather  small  value  to  a  value  of  order  unity.  Since  the  quantities  r\^^  are  the  only  non  vanishing 
elements  of  the  transition  matrix  of  a  sphere,  Eq.  (14)  permits  to  associate  each  resonance  to  one 
and  only  one  of  the  elements  of  the  matrix  S.  For  a  nonspherical  scatterer  the  occurrence  of  a 
resonance  is  no  longer  governed  by  the  behavior  of  the  diagonal  elements  of  the  transition  matrix, 
because  the  off  diagonal  elements  of  S  may  play  an  important  role.  As  a  result,  it  is  rather  difficult 
to  associate  the  observed  resonances  of  a  non  spherical  scatterer  to  a  particular  element  of  the 
transition  matrix  without  a  close  examination  of  the  behavior  of  all  the  relevant  elements  of  S  as 
a  function  of  the  wavelength.  On  the  basis  of  the  preceding  considerations  we  decided  to  label  the 
resonances  both  of  spherical  and  of  non  spherical  scatterers  as  magnetic  and  electric  resonances 
according  to  whether  they  are  associated  to  the  magnetic  part,  or  to  the  electric  part, 

5  appropriate  transition  matrix;  in  this  respect  let  us  remark  that  for  the  scatterers 

that  we  are  going  to  describe  never  occurred  that  a  resonance  should  be  associated  with  the  mixed 

elements  or 

eiemems  or 


A.  Single  spheres 

Our  first  step  has  been  the  calculation  of  the  dependent  scattering  from  the  actual  sphere  and 
its  image  as  a  function  of  the  distance  d  of  the  center  of  the  actual  sphere  from  the  surface.  By 
assuming  that  the  homogeneous  medium  that  fills  the  accessible  half-space  is  the  vacuum  {n  =  1) 
we  found  that,  for  a  homogeneous  spherical  particle  of  radius  p  and  (real)  refractive  index  no  ==  3, 
the  multiple  scattering  processes  between  the  particle  and  its  image  become  negligibly  small  within 
the  whole  range  of  size  parameter  in  which  we  are  interested  when  d  =  lOp.  Thus,  if  the  sphere  is 
at  least  so  far  from  the  interface  we  can  use  the  formulas  that  we  reported  at  the  end  of  Section  3. 
In  particular,  for  the  quantity  Eq.  (24),  a  straightforward  calculation  yields  the  result 

^  (_).-i2^(2/ + 1) + (25) 

m 

on  account  of  Eqs.  (10)  and  (15)  and  of  the  orthogonality  properties  of  the  W  amplitudes. 

We  report  in  Figs.  1  (a),  (b),  (c)  and  (d)  the  plots  of  for  /  <  4,  t;  =  1,  2  and  p  =  1,2, 
as  a  function  of  'dj.  We  notice  in  Figs.  1  (a)  and  (b)  that  at  =  0°  all  the  vanish  for  even 
/  whereas  all  the  do  vanish  for  odd  /.  This  result  was  expected  because  when  dj  =  the 
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only  non- vanishing  amplitudes  are  those  with  in  —  ±1.  Therefore  the  above  result  is  a 

consequence  of  the  fact  that  the  factor  within  square  brackets  in  Eq.  (10)  reduces  to  1  -|- 
Although  the  preceding  argument  is  based  on  the  vanishing  of  the  individual  amplitudes 
at  ??/  =  0°  the  usefulness  of  the  quantity  remains  unaffected.  In  fact,  the  vanishing  of 
for  /  —  2  at  1^/  =  45®,  Fig.  2  (b),  is  due  to  the  sum  over  m  in  Eq.  (24).  The  plots  in  Figs.  1  (c) 
and  (d)  show  that  the  behavior  of  similar  to  that  of  U^\  We  remark  that  again  all  the 

vanish  at  ??/  =  0°  for  even  /  whereas  all  the  vanish  at  the  same  incidence  for  odd  /  and, 
in  particular,  turns  out  to  be  identically  zero:  these  features  were  expected  on  the  ground  of 
the  structure  of  Eq.  (25). 

We  report  in  Figs.  2  (a),  (b)  and  (c)  the  quantity 

=  2^'Im[(-)^~V^;^(ki?,k/)], 

for  r;  =  1,  for  the  sphere  referred  to  above,  in  the  presence  of  the  reflecting  surface,  as  a  function  of 
the  size  parameter  x  —  nkp.  The  respective  values  of  angle  of  incidence  are  di  =  0®,45®,70®,  the 
latter  choice  being  suggested  by  the  fact  that,  according  to  Figs.  1  (a)  and  (b),  at  this  incidence 
none  of  the  quantities  u[]^  vanish  although  several  of  them  assume  a  small  value.  The  quantity 
7i  is  as  meaningful  as  (Text  i  because  it  gives  the  extinction  coefficient  of  a  low  density  dispersion 
of  identical  scatterers^^ .  In  each  of  Figs.  2  we  also  report  the  plot  of  the  quantity 

7  =  2Mm[^r;(k/,k/)], 

for  the  same  sphere  in  the  absence  of  any  substrate.  All  the  resonances  were  classified  with  the 
help  of  the  well  known  formulas. We  stress  that,  in  spite  of  the  indices  of  /,  the  quantity  7  is 
independent  of  the  polarization,  and  that,  for  our  purposes,  7^^  is  quite  comparable  to  7  because 
both  quantities  refer  to  forward  scattering,  according  to  the  discussion  in  Section  3. 

The  strict  correspondence  between  the  resonances  that  disappear  and  the  vanishing  of  the  re¬ 
spective  is  so  evident  that,  in  our  opinion,  no  further  comment  would  be  necessary.  However, 
we  call  the  attention  of  the  reader  on  the  simultaneous  disappearence,  in  Fig.  2  (a),  of  the  two 
resonances  at  a:  =  1.3118,  that  is  associated  to  p  ==  2  and  /  —  1,  and  at  a:  =  1.437,  that  is  associated 
to  p  ==  1  and  1  —  2:  both  peaks  belong,  in  fact,  to  an  odd  value  of  p  +  /.  The  same  mechanism 
explains  also  the  simultaneous  disappearence  of  the  peaks  at  at  x  =  2.2  and  at  x  —  2.28.  The 
former  peak  is  associated  with  p  —  2  and  /  =  3  whereas  for  the  latter  peak  p  =  1  and  /  =  4:  again, 
both  peaks  belong  to  an  odd  value  of  p  -|-  /.  Even  the  resonance  spectrum  for  r}  —  2  strictly  follows 
the  behavior  of  the  corresponding  so  that  we  resolved  not  to  report  the  specific  plot  that,  in 
spite  of  its  significance,  do  not  add  any  further  information  worth  of  a  separate  comment. 

B.  Binary  clusters 

The  lack  of  a  general  theory  applicable  to  non-spherical  particles  makes  rather  difficult  an  unam¬ 
biguous  classification  of  their  resonances.  According  to  Eq.  (14),  the  lack  of  diagonality  of  the 
transition  matrix  prevents  a  meaningful  definition  of  a  function  analogous  to  the  quantity 
that  we  defined  above  for  single  spheres.  Even  in  the  case  of  aggregated  spheres  the  transition 
matrix  is  not  a  diagonal  matrix  so  that  there  is  no  one-to-one  association  of  the  multipole  ampli¬ 
tudes  of  the  exciting  field  to  those  of  the  scattered  field;  nor  there  is  a  simple  relation  between  the 
resonances  of  the  component  spheres  and  those  of  the  aggregate  as  a  whole.  Nevertheless,  Eq.  (15) 
does  not  depend  on  the  shape  of  the  particles  so  that  the  vanishing  of  any  of  the  amplitudes  of 
the  exciting  field  is  expected  to  affect  the  resonances  even  of  a  non  spherical  particle.  To  show 
that  this  is,  indeed,  the  case  we  resolved  to  investigate  the  resonance  spectrum  of  an  aggregate  of 
two  identical  mutually  contacting  spheres  of  radius  p  both  in  the  absence  and  in  the  presence  of  a 
perfectly  reflecting  surface.  In  this  respect  we  recall  that  the  procedure  for  the  calculation  of  the 
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transition  matrix  of  aggregated  spheres  is  outlined  in  Ref.  15.  Even  in  this  case  the  medium  that 
fills  the  accessible  half-space  was  chosen  to  be  the  vacuum  (77.  =  1);  for  the  refractive  index  of  the 
component  spheres  was  assumed  the  unusually  high  value  of  no  =  31.4.  According  to  Newton^^ 
and  to  our  prevous  experience^,  this  choice  makes  the  resonances  both  of  the  component  spheres 
and  of  the  aggregate  as  a  whole  to  occur  at  so  small  values  of  x  =  nkp  that  fully  convergent  values 
of  the  scattered  field  are  obtained  from  Eq.  (14)  for  /  =  /'  =  1  only.  As  a  result,  the  resonances 
of  the  aggregate  as  a  whole  can  surely  be  associated  to  the  multipole  amplitudes  with  /  =  1.  In 
view  of  the  anisotropy  of  any  aggregate  of  spheres  its  orientation  with  respect  to  the  incident  field 
must  be  stated  even  when  no  reflecting  surface  is  present.  The  axis  of  the  binary  aggregates  that 
we  consider  was  chosen  to  be  parallel  to  the  x  axis  of  the  frame  of  reference  that  we  introduced  in 
Section  2.  The  plane  of  incidence,  in  turn,  always  coincide  with  xz  plane.  Even  for  an  aggregate 
of  spherical  scatterers  of  given  orientation  it  is  meaningful  to  introduce  the  quantity  7^  in  strict 
analogy  to  the  definition  that  we  gave  for  spheres  in  Subsection  A. 

In  Fig.  3  (a)  we  report,  for  both  values  of  ?/,  the  quantity  jrj  for  the  cluster  referred  to  above 
in  the  absence  of  any  reflecting  surface  whereas  in  Fig.  3  (b)  we  report  jr)  for  the  case  in  which 
a  reflecting  surface  is  present  and  the  axis  of  the  aggregate  lies  at  a  distance  d  —  lOp  from  the 
surface.  As  in  the  case  of  the  single  spheres,  indeed,  our  calculations  show  that  at  this  distance 
the  multiple  scattering  processes  between  the  actual  aggregate  and  its  image  are  quite  negligible. 
All  the  plots  are  reported  as  a  function  of  x  =  nkp  and  the  angle  of  incidence  is  dj  =  70°. 

Figures  4  (a)  and  (b)  are  identical  to  Figs.  3  (a)  and  (b),  respectively,  except  that  the  angle  of 
incidence  is  7?/  —  0°. 

On  the  whole,  Figs.  3  and  4  present  three  peaks  at  xi  =  0.09629,  X2  =  0.09813  and  xs  =  0.10183. 
Since  for  /  =  1  the  component  spheres  show  a  single  resonance  at  x  =  0.1,  the  plots  in  Figs.  3 
and  4  confirm  our  previous  statement  in  Section  1  that  the  multiple  scattering  processes  within  an 
aggregate  produce  resonances  whose  location  cannot  be  related  to  the  locations  of  the  resonances 
of  the  component  spheres. 

The  dependence  of  the  resonanace  spectrum  of  the  aggregate  on  the  polarization  of  the  incident 
light  can  be  understood  only  through  a  close  examinantion  of  the  features  of  the  transition  matrix 
and  of  the  amplitudes  of  the  exciting  field.  In  this  respect  it  is  important  to  recall  that  the  matrix 
S  does  not  depend  on  the  polarization  so  that  the  dependence  on  the  polarization  of  the  spectra 
in  Figs.  3  and  4  are  entirely  due  to  the  properties  of  the  incident  amplitudes  • 

Now,  at  xi  the  largest  elements  of  S  are  ^  ^  magnetic  resonance 

{p  —  1)  is  expected;  at  x 2  the  leading  elements  are  S[^q^I  q  and  ^ 

electric  resonance  {p  =  2)  is  expected;  finally  at  Xs  the  leading  elements  of  S  are  5*110^1,0 

-1  ^  i  1  suggesting  that  this  resonance  is  a  magnetic  one.  Nevertheless  by 

comparing  Figs.  3  and  4  one  sees  that  not  all  the  possible  resonances  actually  occur.  This  is  due 
to  the  dependence  of  the  amplitudes  both  on  the  polarization  and  on  the  angle  of  incidence 

As  an  example  let  us  discuss  the  behavior  of  the  resonance  at  X2  that  appears  in  Fig.  3  (a) 
for  any  choice  of  the  polarization.  According  to  Eq.  (14)  the  implied  amplitudes  of  the  incident 
field  are  =  (-)” ^  o  =  0-  Therefore,  when  rj  =  1 

the  peak  at  X2  belongs  to  m  =  0  whereas  for  7?  =  2  it  belongs  to  m  =  ±1.  When  the  reflecting 
surface  is  present,  the  implied  amplitudes  of  the  exciting  field  are  ±1  =  0,  q  ^  0  and 

W^2  1  0  “  enough  to  explain  why  the  peak  at  X2  may  appear  in  Fig.  3  (b)  only  for  7;  =  1. 

Wnen  =  0°  the  appropriate  resonance  spectra  are  those  in  Fig.  4.  The  behavior  of  the  peak 

(2)  (2') 

at  X2  is  easily  understood  when  one  considers  that  at  dj  =  0°  we  have  ~  1,1,0  — 

Therefore  a  resonance  can  occur  in  Fig.  4  (a)  for  77  —  2  only  and  cannot  appear  at  all  in  Fig.  4  (b). 

The  behavior  of  the  resonances  at  xi  and  at  X3  that  appear  in  Figs.  3  and  4  can  be  easily 
understood  through  a  quite  similar  analysis. 
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5.  Conclusions 


The  results  that  we  presented  in  Section  4  suggest  that  the  presence  of  a  perfectly  reflecting  surface 
may  be  a  useful  tool  for  the  interpretation  of  the  resonance  spectra  from  scattering  particles. 
However,  the  behavior  of  the  spectra  from  spherical  scatterers  must  be  carefully  distinguished 
from  the  behavior  of  the  spectra  from  particles  of  more  general  shape. In  fact,  the  diagonality 
of  the  transition  matrix  of  spheres  simplifies  the  association  of  the  observed  resonances  to  the 
multipole  amplitudes  of  the  exciting  field:  thus  each  resonance  can  easily  be  classified  as  electric 
or  magnetic  and  attributed  to  the  appropriate  value  of  /.  Further  insight  into  the  behavior  of  the 
resonances  as  a  function  of  direction  of  propagation  and  of  the  polarization  of  the  incident  wave  is 
supplied  by  the  quantity  that,  as  we  stressed  in  Section  4,  can  be  defined  only  for  spheres. 

The  interpretation  of  the  resonances  of  particles  of  general  shape  is,  on  the  contrary,  rather 
difficult.  In  fact,  since  the  transition  matrix  of  such  particles  is  non-diagonal,  several  elements  of 
S  csin' contribute  to  the  same  resonance.  As  a  result,  it  is  not  easy  to  establish  a  correspondence 
between  the  observed  resonances  and  the  multipole  amplitudes  of  the  exciting  field. 

Even  when  the  particles  of  interest  are  or  can  be  modeled  as  clusters  of  spherical  scatterers  there 
is,  in  general,  no  evident  relation  between  the  resonances  of  the  aggregate  as  a  whole  and  those  of 
the  component  spheres.  In  fact,  the  discussion  in  Section  4  on  the  resonance  spectrum  of  binary 
clusters  has  been  made  possible  only  by  choosing  a  small  value  for  the  radius  and  a  rather  high 
value  for  the  refractive  index  of  the  component  spheres.  Due  to  our  choice  the  resonances  of  the 
aggregate  occur  for  x  <C  1  so  that  the  convergence  of  the  scattered  field  requires  to  include  into 
the  multipole  expansions  terms  with  /  =  1  only.  For  realistic  values  of  the  refractive  index  the 
resonances  of  the  cluster  occur  for  higher  values  of  x  so  that  convergence  of  the  scattered  field 
requires  to  include  into  the  multipole  expansion,  Eq.  (13),  terms  with  higher  values  of  /:  this, 
in  turn,  allows  further  resonances  associated  with  higher  values  of  /  to  appear.  In  spite  of  these 
difficulties,  the  simple  example  that  we  reported  in  Section  4  suggests  that  a  careful  analysis  as 
a  function  of  ??/  and  of  the  polarization  both  in  the  presence  and  in  the  absence  of  a  reflecting 
surface  may  give  useful  information  for  the  interpretation  of  the  resonance  spectra  of  non-spherical 
particles. 
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Fig.  1.  Plot  of  the  quantity  for  I  <  4.  In  (a)  /?  —  1  and  r}  =  1;  in  (b)  p  =  2  and  77  =  1;  in  (c) 

p  =  1  and  7;  =  2;  in  (d)  p  =  2  and  77  =  2.  In  this  paper  we  adopted  the  convention  that  p  =  1,  2  classifies 
the  multipoles  as  magnetic  and  electric,  respectively;  in  turn  77  =  1,2  indicates  that  the  polarization  is 
parallel  or  orthogonal  to  the  plane  of  incidence,  respectively. 


Fig.  2.  7t7  (solid  curve)  and  7  (dotted  curve)  for  a  homogeneous  sphere  of  radius  p  and  refractive  index 

770  =  3  as  a  function  of  x  =  nkp  for  r]  =  1.  The  medium  that  fills  the  accessible  half-space  ia  assumed  to  be 
the  vacuum  (n  =  1).  The  distance  between  the  center  of  the  sphere  and  the  reflecting  surface  is  d  —  lOp. 
The  angle  of  incidence  is  -dj  =  0°  in  (a),  di  =  45®  in  (b)  and  di  =  70®  in  (c).  The  resonances  are  classified 
according  to  the  scheme  (p,0«  )  where  n  distinguishes  different  resonances  with  the  same  value  of  p  and  1. 


Fig.  3.  7t7  for  the  aggregate  of  two  identical  mutually  contacting  spheres  of  radius  p  and  refractive 

index  tiq  =  31.4  as  a  function  of  x  =  nkp.  The  medium  that  fills  the  accessible  half-space  ia  assumed 
to  be  the  vacuum  (n  =  1).  The  axis  of  the  aggregate  lies  in  the  xz  plane  and  is  parallel  to  the  x  axis. 
The  plane  of  incidence  coincides  with  the  xz  plane.  The  angle  of  incidence  is  di  =  70®.  The  solid  and 
the  dotted  curves  refer  to  polarization  parallel  and  orthogonal  to  the  plane  of  incidence,  respectively.  In 
(a)  no  reflecting  surface  is  present,  whereas  in  (b)  a  reflecting  surface  that  coincide  with  the  xy  plane  is 
present  and  the  distance  between  the  axis  of  the  aggregate  and  the  surface  is  d  =  lOp.  We  notice  that  the 
spike  at  X  =  0.09813  in  Fig.  3  (a)  appears  for  both  choices  of  the  polarization,  although  this  is  not  easily 
discernible  on  the  scale  of  the  figure. 


Fig.  4.  Same  as  Figure  3  except  that  the  angle  of  incidence  is  now  d/  =  0®. 
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